Phase Transitions of Soluble Surfactants at a Liquid- Vapor Interface 
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Although medium chain length insoluble amphiphiles are well known to form gaseous and liquid 
expanded phases on an air/water interface, the situation for the soluble case is less clear. We 
perform molecular dynamics simulations of model surfactant molecules dissolved in a bulk liquid 
solvent in coexistence with its vapor. Our results indicate a transition in both soluble and insoluble 
surfactants: a plateau in surface tension vs. surface coverage, whose instantaneous configurations 
display two phase coexistence, along with correlation functions indicating a transition to gaseous to 
liquid-like behavior. 

PACS numbers: 79.20.Rf, 64.60. Ht, 68.35.Rh 



Surfactants are amphiphilic molecules which contain a 
polar head group and a nonpolar tail, which usually con- 
sists of a chain of hydrocarbon groups. At an air/water 
interface surfactants arrange themselves with their polar 
group immersed in water, interacting by dipolar forces, 
while the hydrocarbon tails are displaced outward into 
the air. The presence of surfactants at the air/water in- 
terface lowers the free energy, and in turn the surface ten- 
sion, usually as a non-increasing function of the concen- 
tration of surfactant in the liquid The degree of ag- 
gregation and molecular ordering of surfactants adsorbed 
at an air/solvent interface is reflected in the appearance 
of distinct surface phases as a function of surfactant con- 
centration. In addition to the obvious two dimensional 
Gas (G) and disordered liquid phases (more precisely, the 
Liquid Expanded LE phase), it is possible to find con- 
densed mesophases and semi-solid crystalline-like phases 
in monolayers ||||. 

This paper focuses on the determination of soluble sur- 
factant adsorption isotherms simultaneously with their 
phase behavior and surface tension variation, and in par- 
ticular the G/LE transition, using Molecular Dynamics 
(MD) simulations. In the context of insoluble surfac- 
tants, this approach has been used extensively to study 
the structure of monolayers restricted to the surface of a 
substrate obtaining information both complemen- 

tary to and in semi-quantitative agreement with exper- 
iment on Langmuir films. For soluble surfactants there 
is emerging experimental evidence for the existence of 
gaseous, liquid, and condensed surface phase transitions 
, but no accompanying theory or simulation. There are 
MD simulations of soluble surfactant monolayers S , but 
these do not study the surface-bulk equilibrium which 
gives rise to the surface phases, or the transitions be- 
tween them. In contrast, MD simulations in liquid- liquid 
systems ^] have studied the equilibrium between the liq- 
uid phase and monomeric and surfactant aggregates in 
the bulk, but not surface phase transitions. 



Our simulations model surfactant molecules placed in 
a solvent in equilibrium with its vapor, which then mi- 
grate to the surface to a degree controlled by the choice 
of interaction potential. We observe that as their ini- 
tial concentration increases, the surfactant structure on 
the liquid/ vapor interface ranges from a gas phase con- 
taining isolated molecules or small clusters, to a coexist- 
ing mixture of clusters of various sizes, to a single disor- 
dered spanning liquid cluster. Examination of positional 
and tail-orientation distribution functions show a gradual 
transition to liquid-like ordering, and the isotherm of sur- 
face tension vs. surface concentration exhibits a plateau 
indicative of a phase transition. We have studied both in- 
soluble and soluble surfactant systems, obtaining rather 
similar behavior, but in this Letter we concentrate on the 
latter. Further details will be presented elsewhere flQfl . 

Our calculations are based on standard MD tech- 
niques using Lennard- Jones (LJ) interactions, in an NVT 
ensemble jll]]. The potential between any two atoms 
of type i and j separated by distance r is (r) = 
Ae[(a/r) 12 — dj (a/r) 6 ]. Here e is an energy scale, a 
is approximately an atomic diameter, and the charac- 
teristic microscopic time unit is t = y/ (ma 2 /e), with 
atomic mass m. The adjustable coefficients CV, de- 
termine the interactions between the various molecular 
species and control the chemistry of the system. The 
surfactant atoms are bound in chain molecules using the 
Finitely Extensible Nonlinear Elastic (FENE) potential 
H V FENE = -\Kr% In [1 - (r/r ) 2 ] which acts only 
between adjoining atoms on a chain; ro is the maximum 
bond length and if is a spring constant. The system is 
maintained at a temperature T — 0.9e/k by a constant 
kinetic energy thermostat. We have worked with a sys- 
tem of 11520 atoms in total, with surfactant chains of 
length 6. The atoms move in a three dimensional box 
of size 20.5 ax 20.5 a x 68.4 a, with periodic boundary 
conditions in all directions, and form a horizontal liquid 
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slab with vapor above and below, with two (statistically) 
planar liquid/vapor interfaces. A snapshot of a typical 
configuration is shown in Fig. [l]a. 




FIG. 1. Snapshots of the simulated system; (a) all atoms 
in a soluble case, (b) surfactant molecules only in an insoluble 
case, (c) surfactant molecules only in a soluble case. 

The surfactant molecules are initially distributed at 
random in the liquid, and then migrate to the liquid- 
vapor interface to a degree determined by the choice of 
Cij , so that the solubility is indirectly under control. An 
insoluble surfactant is achieved by choosing a strong at- 
traction between head group and solvent, an antipathy 
between surfactant tail and solvent, and a strong cohe- 
sion among the molecules in the solvent. The full set of 
interaction coefficients is shown in Table I, where the sub- 
script convention is 1: hydrophilic head, 2: hydrophobic 
tail and 3: solvent. A typical snapshot of the insoluble 
surfactant phase alone is shown in Fig. [l]b. For a soluble 
surfactant system, we modify the interactions to allow 
surfactant to remain in solution. A first ingredient is to 
reduce the antipathy between the surfactant tails and the 
solvent, by increasing C23, and the second is to reduce 
the cohesion of the solvent, by reducing C 33 , and the full 
set is shown in Table II, 

The surface tension is computed from the standard 
microscopic expression for a planar interface Jt3| , 7 = 

\ (J2i<j Tzj 2r* ij dV dr z ] j) \ where the sum runs over all 
interacting pairs of atoms, V is the full interatomic po- 
tential, A is area of the interface, and the angle brack- 
ets denote a time average. The insoluble surface tension 
isotherm is given in Fig. ^. There is a clear reduction 
for small amounts of surfactant, followed by an approxi- 
mate plateau until a critical concentration is reached, and 
then a rapid decrease. (At still higher surfactant concen- 
tration we observe buckling of the interface or micelle 
formation.) The structure of the surfactant layer may 



be inferred from the pair distribution function g(r) p4[ , 
here taken to be the probability that two head groups on 
a given interface are separated 



Table I: Insoluble. 



Table II: Soluble. 
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by a distance r. As we see in Fig. [}]a, in the very low 
coverage region g(r) is a flat line typical of a gas, but as 
the coverage increases up to the end of the plateau, the 
curve changes to that of a typical liquid, with a promi- 
nent nearest-neighbor first peak, a rather broad second 
peak exemplifying short range order, and some further 
structure at larger r as well. This simulation therefore 
reproduces the well-documented G/LE phase transition 
for insoluble surfactant. 
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FIG. 3. Pair distribution function for surfactant monolay- 
ers on the liquid-vapor interface: (a) insoluble, (b) soluble. 

For the soluble surfactant system, a representative 
snapshot of the positions of surfactant molecules is in 
Fig. Qcj. The measured isotherm is shown in Fig. ^|, and 
has the same qualitative features as in the insoluble case. 
In this figure, the surfactant surface coverage is deter- 
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mined by considering the component density profiles as 
a function of distance normal to the interface, Fig. ||. 
The upper curve is the solvent profile, while the lower 
one refers to the tail groups, and the surface coverage T 
is taken to be the total area under the bump to the left 
of Z\. 
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FIG. 4. Interface profile obtained in the simulation, and 
the Gibbs construction. 

To obtain the surfactant adsorption isotherm in a ther- 
modynamically unambiguous manner, we use the density 
profiles to extract Gibbs' "surface excess concentration" 
T e We first locate a Gibbs dividing surface, z , which 
serves as a convenient definition of the solvent interface. 
Referring to Fig. ||, the location of this surface is defined 
by the condition that the areas of the shaded regions A 
and B are equal, thereby requiring the surface excess of 
solvent to vanish. The bulk surfactant concentration is 
defined by extrapolating the average bulk density S in 
the bulk up to the Gibbs dividing surfaces, 
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FIG. 5. Surface concentration vs. Bulk concentration. 

and the surface excess concentration of surfactant F e 
is then the remainder, - the area of the hatched region 
in the figure. (We assume that there is no surfactant in 
the vapor phase.) Repeating this procedure at each value 
of bulk concentration, we obtain the variation of surface 
excess coverage with total concentration shown in Fig. |^: 
a monotonic increase, with a sigmoidal shape. Note that 
the steeply rising part of the curve occurs precisely in the 



plateau region of the isotherm, which we will argue cor- 
responds to the G/LE phase transition. In the thermo- 
dynamic limit one would expect a vertical line there, as 
surfactant is expelled from solution to form larger liquid- 
like clusters on the interface, but we observe the smooth- 
ing and finite slope characteristic of phase transitions in 
finite-sized systems. 

Returning to the soluble isotherm, note that the crit- 
ical densities of the gas and liquid phases for the G/LE 
transition (the end points of the plateau) occur near the 
same values of T found for the insoluble case, presum- 
ably because this value corresponds to the onset of clus- 
tering, and is determined largely by the attraction be- 
tween surfactant molecules, which was not changed. The 
pair distribution function for systems at the start, mid- 
dle and end of the surface tension plateau are shown in 
Fig. ||b, and generally resemble the insoluble case. The 
characterization of these phases is facilitated by examin- 
ing snapshots of the molecular configuration, and three 
typical examples corresponding to the three surface cov- 
erage values where g(r) was plotted are shown in Fig. ||. 




FIG. 6. Snapshots of soluble surfactant molecules at an in- 
terface: (a) gas,(b) gas-liquid and (c) liquid expanded phases. 

At low surface concentration, Fig. ||a at T = 0.014, one 
sees a gaseous state with isolated molecules and a weak 
tendency to bind into pairs. At higher coverage T = 0.062 
in the middle of the plateau, Fig. ||b, we see a spanning 
cluster of surfactant molecules, accompanied by isolated 
molecules. Finally at the end of the plateau, T = 0.097 
in Fig. ||c, there is a single surfactant monolayer. 

Further evidence for phase transitional behavior in the 
plateau region is obtained by studying the degree of or- 
dering of the surfactant tails. A simple quantitative mea- 
surement is the histogram of the tilt angle between the 
surfactant tail and the normal to the interface, shown in 
Fig. [?]. The orientation of a molecule is defined as the 
direction of the eigenvector of its moment of inertia ten- 
sor with the smallest eigenvalue ||. We see that as the 
surface coverage increases across the plateau, molecules 
are increasingly likely to be oriented normal to the inter- 
face. An even more dramatic variation is seen in the his- 
togram of end-to-end lengths, shown in the inset to Fig. |t] 
- as coverage increases, the molecules are increasingly 
extended as they pack together. Note that because the 
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molecular interaction corresponds to freely-jointed chains 
with no bond-bending forces, the surfactant molecules 
need not form straight rods, and we do not observe com- 
plete correlation of the tails seen in MD simulations of 
condensed phases in Langmuir monolayers 0. Further- 
more, since the tails are now partially soluble they are 
more likely to have solvent around them, which tends to 
decorrelate their orientation. 
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FIG. 7. Probability distributions of the tilt angles, at vari- 
ous values of surface concentration. The inset shows the corre- 
sponding histogram of the length of the molecules' end-to-end 
vectors. 
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Lastly, we can go further in this vein and define a pair 
angular correlation function: the probability distribution 
for the mean relative angle between two molecular tails 
as a function of the head group separation r. Results are 
shown in Fig. ^: at low coverage, there is little corre- 
lation between the tails, but the likelihood of alignment 
increases with T, as the molecules pack together and the 
tail are attracted to each other. 

In summary, we have presented MD simulations of 
soluble amphiphilic systems, containing simple chain 
molecules in a monatomic bulk solvent, which assem- 



ble from solution into a surfactant layer at the liquid- 
vapor interface. As a function of surface coverage, these 
systems exhibit a G/LE phase transition which can be 
characterized by plateaus in the surface tension isotherm, 
visible changes in the molecular configurations, and cor- 
responding transitions in various correlation functions. 
The behavior of insoluble and soluble systems is found 
to be very similar, suggesting that at least this surface 
ordering phenomenon is not highly sensitive to the dy- 
namics in the bulk solution, but rather an effect of surface 
packing. 
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